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ABSTRACT 

Context. The bulk of solar flare emission originates from very compact sources located in the lower solar atmosphere and seen in 

various wavelength ranges: near optical, UV, EUV, soft and hard X-rays, and gamma-ray emission, yet very few spatially resolved 

imaging observations to determine the structure of these compact regions exist. 

Aims. We investigate the above-the-photosphere heights of hard X-ray (HXR), EUV and white-light (6173 A) continuum sources in 

the low atmosphere and the corresponding densities at these heights. Considering collisional transport of solar energetic electrons we 

also determine where and how much energy is deposited and compare these values with the emissions observed in HXR, EUV and 

continuum. 

Methods. Simultaneous EUV/continuum images from AIA/HMI on-board SDO and HXR RHESSI images are compared to study a 

well observed gamma-ray limb flare. Using RHESSI X-ray visibilities we determine the height of the HXR sources as a function of 

energy above the photosphere. Co-aligning AIA/SDO and HMI/SDO images with RHESSI we infer, for the first time, the heights and 

characteristic densities of HXR, EUV and continuum (white-light) sources in the flaring footpoint of the loop. 

Results. 35-100 keV HXR sources are found at heights between 1.7 and 0.8 Mm above the photosphere, below the 6173 A continuum 

emission which appears at heights 1.5-3 Mm, and the peak of EUV emission originating near 3 Mm. 

Conclusions. The EUV emission locations are consistent with energy deposition from low energy electrons of ~ 12 keV occurring in 

the top layers of the fully ionized chromosphere/low corona and not by > 20 keV electrons that produce HXR footpoints in the lower 

neutral chromosphere. The maximum of white-light continuum emission appears between the HXR and EUV emission, presumably 

in the transition between ionized and neutral atmospheres suggesting free-bound and free-free continuum emission. We note that the 

energy deposited by low energy electrons is sufficient to explain the energetics of optical and UV emissions. 
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1. Introduction 

Solar flare emission is now observed virtually over the whole 
electromagnetic spectrum from radio frequencies as low as 10 - 
100 MHz up to photon energies as high as a few hundred MeV. 
The bulk of this emission originates from the chromosphere, a 
small part of the solar atmosphere which is narrow in height and 
much smaller than the size of a flaring region. Being only about 
< 3 Mm (< 4'') across, the chromosphere presents a formidable 
observational challenge in both angular and temporal resolution 
especially for transient phenomena like solar flares. During so- 
lar flares large numbers of energetic electrons are accelerated 
in the solar corona and eflfectively release their energy in the 
dense chromosphere. These energetic particles are believed to 
be responsible (although not always directly) for HXR, EUV, 
continuum "white-light" (WL), infrared and radio emissions. 
On the basis of temporal correlation, th e interconnection be- 
tween these emissions was early realized ( Najita & Qrrall|1970t 
'Svestka'l970). Due to the lack of height-resolving observations 
spectroscopic data have been m ostly used to id entify the struc- 
ture of the emitting region (e.g. Canfield| 1974| ). Therefore, the 
actual mechanisms of WL emission and the relation to HXR 
emitting electrons is still subject of debate ([Neidig 1989; Neidig 
et al.|1993t[Sylwester & Sylwester|20QQl [Matthews et al.|2003[ 



[Potts et al.|20T0l|Watanabe et al.|2010t|Kretzschmar|2011| ). It is 
not even clear whether the emission is optically thick or thin or 
both. Thus the characteristic heights of continuum WL emission 
in various models are placed from the photosphere to the up- 
per chro mosphere. Ob servations supporting both (Fang & Ding 
1995; X u et al.[[2006|), a photospheric origin (e.g. Boyer et al 



|1985f |Ding et al.[[ 19991 [Chen & Ding|[2QQ6l ) presumably due 



to radiative back-warming, and generation of WL emission in 
the higher temperature regions of the upper chromosphere (e.g. 
.Machado & Rust, 1974^, Hudson, 1 972 j have been made. The ma- 
jor observational challenge here is to have simultaneous HXR, 
UV and WL images with sufficiently high angular resolution. 
The relative horizontal positions of HXR and WL flares have 
been investig ated in detail and, although generally been found 
to agree (e g.|Krucker et al.|2011[ [Fletcher et al.|2007[ [Hudson 



et al.[[2006| [Metcalf et al.[[2b03[ ), the uncertainty in the relative 
position was large due to limited absolute pointing accuracy of 
partial solar disk observations such as TRACE. Thus it was not 
possible to quantify the relative height of the sources above the 
photosphere in previous studies. With RHESSI ( Lin et al.|2002 l 
data, the height structure of HXR sou rces has become accessi- 
ble for detailed observational studies (Aschwanden et al. "20021 
Kontar et aL[[2008, Kontar & Jeffrey. 2010, Battaglia & Kontar] 
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201 It [Saint-Hilaire et al.[[2010[ ). These observations demon- 
strated that the higher energy HXR footpoints originate from the 
chromosphere from progressively lower heights, in good agree- 
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Fig. 1. Top: RHESSI lightcurves of the main flare phase in dif- 
ferent energy bands: 6-12 keV (green), 12-25 keV (red), 25-50 
keV (blue), 50-100 keV (purple). Bottom: Normalized AIA and 
HMI emission from location of the southern HXR footpoint as 
function of time in several wavelengths (see legend). 



ment with coUisional downward propagation of energetic elec- 
trons (see Brown et"aLl|2002 for details). The newly available 
SDO conti nuum images from HMI ( Wachter et al.|[2QTl ) and 
AIA EUV ( [Lemen et al.|2011 ) observations allow to scrutinize 
the spatial structure and more importantly the height structure 
of various emissions associated with HXR footpoints. In this 
paper we present the first spatially resolved simultaneous ob- 
servations of HXR, EUV and WL emission from a footpoint of 
a solar limb flare observed with RHESSI, SDO/AIA and HMI. 
Using the recently developed hard X-ray visibility analysis tech- 
nique ( Hurford et al..2002l [Schmahl et"ari[2007) we determine 
the characteristic heights for the diff'erent types of emission us- 
ing the method presented by |B attaglia & Koiit^ ( |20 1 1 1 ) ; [Kontar 
|etal.| ([2010, 2008). The results are consistent with the EUV and 
continuum WL being produced above the HXR footpoints by 
lower energy electrons (-12 keV) and hence suggest not a pho- 
tospheric but upper chromospheric origin of WL emission. 



2. Height and density measurements in HXR and 
EUV 

The GOES M3.5 limb flare happened on February 24th 2011 
with three well-pronounced main HXR (above ~ 30 keV) peaks 
between 07:29 and 07:33 UT and an associated filament eruption 
to the south of the flaring site. Figure [T] (top) shows the RHESSI 
lightcurves in different energy bands of the thermal and non- 
thermal emission. The attenuator state was 1 during the main 
phase of the flare and the live-time was better than 90% dur- 
ing the whole period of the observations, hence making the ef- 
fects of pulse pile-up negligible. The flare appeared in active re- 
gion ARl 1 163 and is well observed in enhanced EUV emission 
visible in all AIA filters, and also showed 6173 A continuum 
emission observed by HMI. The continuum emission (with a 
maximum enhancement of ~ 15% above the photospheric back- 
ground) shows two footpoints and follows the X-ray time pro- 
file. Unfortunately, AIA images became saturated in most wave- 
lengths before the soft X-ray peak (2nd HXR peak) of the flare. 
Because of this saturation, we focus the analysis on the time of 
the first HXR peak. During the time around 07:30:10 UT, AIA 
was not saturated in 94 A, 193 A, 304 A and 335 A filters near 
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Fig. 2. Top: AIA 171 A image overlaid with RHESSI soft (red) 
and hard X-ray emission (blue). Bottom: AIA diff'erence images 
in 3 wavelengths and HMI image (bottom right). Contours at 
50, 70 and 90 % of the maximum emission in RHESSI CLEAN 
images between 07:30:00 and 07:30:40 UT at 6-12 keV (red) 
and 25-50 keV (blue) are overlaid. 



the location of the southern footpoint. To study the flare related 
enhancement in EUV and continuum emissions the image from 
one minute before the analyzed time was subtracted. Figure |2] 
displays difl'erence images in AIA 94, 193 and 304 A wave- 
lengths and the HMI continuurrQ For the RHESSI analysis a 40 
second time interval at the first peak of the flare from 07:30:00 
UT was chosen (see Fig.[T]). The contours in Fig. [2] indicate the 
RHESSI emission between 07:30:00 and 07:30:4 UT at 50, 70 
and 90 % of the maximum in CLEAN images ( [Hurford et al 



2002). The images show a coronal (SXR) source above the limb 



and two footpoints (HXR) sources (north and south). During the 
first peak the southern footpoint source was much more intense 
than the northern footpoint. Thus we focused the height-analysis 



^ Animated images of the flare are available online: 
|http ://w w w. astro . gla. ac .uk/user s/mbattaglia/20 1 1 0224 _online -material/] 
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on the southern footpoint. We note that in the subsequent peaks 
the northern footpoint reached the same intensity as the southern 
footpoint. The RHESSI spectrum suggests emission well beyond 
200 keV into the y-ray range. During the analyzed time interval, 
the HXR spectrum is well fitted with a thermal component and 
a single thick- target power-law with spe ctral index S = 3.9 . The 
flare was also observed by Eermi/GBM (Meegan et al."2009 ) and 
the corresponding spectrum from EERMI data was found virtu- 
ally identical to the RHESSI spectrum. The time evolution of the 
EUV emission at the position of the southern HXR footpoint can 
be seen in Eig.[T] (bottom). The total flux in a region of 5 by 10 
arcsec extent around the position of the southern HXR footpoint 
is shown. 

Because of the limb location of the flare the height analy- 
sis is possible. However, SDO observations alone cannot infer 
the heights of the sources in difl'erent EUV wavelengths or for 
white-light emission. Here, utilizing RHESSI observations, we 
find the density structure of the chromosphere and importantly 
the reference height (projected radius of the point at t he photo- 
spheric level beneath the HXR sources) as described by Battaglia] 
|&Kontar(20ir ) and |Kontar et al.| ( |25TQ|[2008j . Eollowing the 
method described by those authors we use visibility forward 
fitting to find the position of the southern footpoint as a func- 
tion of HXR energy in several energy bands (30-40 keV, 40- 
55 keV, 55-80 keV, 80-130 keV). The energy bands were cho- 
sen large enough to guarantee good count statistics for reliable 
fits. Both the southern and the northern footpoint are fitted with 
circular Gaussians. This assures that the emission of the north- 
ern footpoint is properly accounted for, even if it is not intense 
enough to result in fit parameters with acceptable errors. The 
flare morphology was such that the radial direction (direction 
of electron propagation along the loop) can reasonably be seen 
as along the x-direction in the images (see Eig. [2]). An expo- 
nential density profile (hydrostatic atmosphere at lower heights), 
n(h) = Ties + nph exp(-(/z - r^e/)/^o). was fitted to the positions 
found by visibility forward fitting, where ^^^ = 4.6 x 10^^ cm"^ 
is the (constant) loo p density, riph = 1.16 x 10^^ cm"^ the pho- 
tospheric density ( [Vernazza et al. [T98T] ) which corresponds to 



TSOOn 



1, 



and h is the height above the photosphere q] The 
loop density rics was assumed to be the same as the coronal 
source density determined from the emission measure given by 
the spectral fit and the size of the coronal source in RHESSI 6- 
12 keV images. Eitting the HXR energy-position relation with 
the above density model we find a scale height of /zq = 21 1 ± 46 
km and a reference height for the radial distance between the 
solar disk center and the point under the footpoints at the pho- 
tosphere of rref = 929.4 ±0.3 arcsec. The photospheric refer- 
ence distance rref is used to find the absolute height of the EUV 
and WL sources above the photosphere. Using the HXR source 
size as the size of the magnetic flux tube, we selected a stripe 
along the x-direction of 5 arcsec width in the y-direction over 
the position of the southern RHESSI footpoint. At each horizon- 
tal distance x, the AIA and HMI counts were summed over the 
y-direction. The height of the emission is simply the x-position 
minus rref as found above. This is shown in Eig. l3] Proper co- 
alignment of the images is crucial for this study. The RHESSI 
disk center is known to better than 0.2 arcsec accuracy. Using 
SDO full disk images the disk center of AIA and HMI can be de- 



^ In this work we use a fit to the full expression for the HXR flux 
maximum and not the simplified version given by Eq. 7 of [Brown et al.| 
( |2002| ). The latter is likely to und erestimate the densities in a sharply 
changing chromos phere. See also [Battaglia & Kontar| ( [20TT] ); |Kontar| 
inL] ( [20T0|[2008] ) 
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Fig. 3. Top: Height of HXR emission at the lowest (red verti- 
cal lines) and highest (blue vertical lines) imaged energies, EUV 
emission at difl'erent wavelengths (see legend) and WL relative 
enhancement where the black line indicates the height of max- 
imum enhancement. Arrows indicate the corresponding uncer- 
tainties. Middle: Energy deposition rate as a function of height 
for two difl'erent cases of electron cutofl' energy. Bottom: Eitted 
density model (thick solid line). The thin solid lines give the 
lower and higher limits of the density due to the uncertainty in 
the fit of the scale height. The model from Vernazza et al.| ( |198l| 
is given for comparison. 



termined. This position of solar-disk center is found to deviate by 
only -0.39/0.67 arcsec in x/y direction that of RHESSI, which, 
although small, were taken into account for the co-alignment. 
The larger uncertainty is related with the roll-angle of SDO im- 
ages (rotation around the disk center). To investigate the eff'ect, 
the HMI image was rotated by 0.1 degree relative to Sun cen- 
ter. This uncertainty translates to mostly j-direction uncertainty 
of < 1.6 arcseconds and < 0.5 for the x-direction for the limb 
event under study. Eor the height measurements of the EUV and 
HMI emissions relative to the RHESSI emission, we have taken 
this uncertainty into account and shown in Eig. [3] We find that 
the RHESSI sources at 35 - 100 keV originate from heights be- 
tween 1.7 and 0.8 Mm. The WL emission peaks at 2.5 Mm, 
which is ~ 1 Mm above the 30 keV HXR emission, suggest- 
ing it originates from the upper chromosphere. Using the den- 
sity model found above and a power-law electron injection rate 
Fo(Eo) = FnormEQ^ [clcctrons s"^ keV"^] found from the spa- 
tially integrated RHESSI X-ray spectrum, we can compute the 
energy deposition rate [ergs s"^ Mm"^] as a function of height 
( [Brown|1973] ), 



dU 
dx 



KlAeHnnih) + A,,^,(/z)] 



X 



F(E,h) 



dE 



(1) 



where F(E, h) is the electron rate spectrum C Brown|1971 ), K = 
27re^ with AeH and Aee the Coulomb logarithms for electron- 
neutral and electron-electron collisions, and nnih) and rieih) are 
the hydrogen and electron number densities respectively. The 
resulting energy deposition rate as a function of height given 
by Equation ([T]) is shown in Eig. [3] where we assumed a step 
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change in ionization from fully ionized to neutral at a height of 
2 Mm. As expected, the deposition rate is dependent on the low 
energy cut-off for the injected electron spectrum Fq(Eq). From 
the thick target fit to the RHESSI spatially integrated HXR spec- 
trum the cut-off* energy is around Ec = 12 keV and the spectral 
index 6=3.9. This leads to a total power of ~ 2.7 x 10^^ ergs"^ 
However, the cut-off* is generally unknown as the low energy 
part of X-ray spectrum is dominated by thermal emission, so 
lower cut-off" values are possible and likely (e.g. Emslie 2003 [ 
[Hannah et al.|2009] ). However, even with Ec = 12 keV, the en- 
ergy deposited appears much larger than emitted in other wave- 
lengths. Eor the southern footpoint, an enhancement of ~ 5% 
over the total solar irradiance of 3.4 x 10^^ergs"^arcsec"^ over 
a source area of ~ 20 arcsec^ leads to a total power emitted 
of 3.4 X 10^^ erg s"^ The energy in the non-thermal electrons is 
therefore up to an order of magnitude larger than what is needed 
to power the WL emission, which often constitutes the bulk of 
ffare emitted energy (Neidig 1989). We note that the WL en- 
ergetics obtained using 6173 A continuum is found to be lower 
than using TRACE wavelengths near the 1700 A continuum, and 
assuming black body emission to deal with the TRACE response 
(e.g. [Fletcher et al" [2007| ), but close to the values reported by 



Neidig ( 1989 ) based on spectroscopy above 2500 A. Indeed, var- 



ious o ptical bands show d iff^erent levels of enhancement during 
flares ( Kretzschmar|20 1 1 1 ) . In addition, the maximum of 6173 A 
continuum is reached later in the flare, after the analyzed time 
interval. 



energy deposition rate. In this case, the maximum energy depo- 
sition of electrons with energy < 12 keV will be near the top of 
the loop. Should the loop density be lower or the acceleration 
region extended, the lower energy electrons would be deposit- 
ing their energy deeper. The other interesting aspect of this flare 
is that HXR emission above 20 keV is well fitted by a single 
power-law. Non-uniform ionization of plasma produces a break 
at the energy corresponding to the column depth of the transi- 
tion region, which is observed in other ff ares ( [Kontar et al. 2002[ 
Su et al.|2011| ). The absence of such break points to the transi- 
tion region above the stopping depth of 20 keV electrons again 
is consistent with the height measurements. The observation of 
strong emission from heights > 2 Mm suggests that a substantial 
part of the WL continuum is formed in an optically thin or finite 
optically thick region with ionized or partially ionized plasma in 
the upper chromosphere favoring free-bound and free-free emis- 
sion of the 6 173 A continuum. This is additionally supported by 
the height of the AIA emission at 94, 335 and 193 A (represent- 
ing temperatures > 10^ K) originating from around 3 Mm and 
above the WL emission peak. Finally, it is likely that the WL 
is powered by lower energy electrons ~ 12 keV, which do not 
penetrate deep into the solar chromosphere. 
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3. Discussion and Conclusions 

Our observations suggest that the WL emission originates from 
the upper chromosphere or lower transition region with a den- 
sity of between lO^^cm"^ and lO^^cm"^. These upper and lower 
limits of the density are due to the uncertainties in the fit of the 
density scale height. Figure [3] (bottom) displays the fitted density 
model, as well as the upper and lower limits for the density. The 
arrows in Fig. |3] (top) illustrate the uncertainty in the RHESSI 
source heights resulting from the uncertainty in the density scale 
height. The peak of WL emission is found about 1 Mm higher 
than the position of the lowest observed RHESSI energy range 
(30 - 40 keV). Note that despite all the uncertainties in the height 
due to the density fit or due to relative pointing error, the WL 
source is still above the 30-40 keV RHESSI sources. Due to the 
presence of the SXR coronal emission in the X-ray spectrum, 
the non- thermal HXR spectrum below ^ 20 keV is not known. 
However, extrapolating the height function found in Section |2] 
to lower energies, a source at 12 keV would be seen at about 2 
Mm. This coincides with the height of the maximum energy de- 
position as found from Eq.[T]for cut-off' energies of 12 keV (Fig. 
l2]). The observed position of the maximum of the WL emission 
would then hint towards energy deposition by the low energy 
part of the non-thermal electron distribution, contrary to the of- 
ten made assumption of electron cut-off' energies in the range 
of 20 keV. This also suggest that WL and HXR sources above 
30 keV should be spatially separated by < 1 Mm with decreas- 
ing separation in events towards the disk center. However, we 
note that the height of the maximum energy deposition, espe- 
cially in the case of low initial cut-off' energies, is inffuenced 
by the loop density and the length of the loop. Assuming that 
the acceleration of the particles happens in the center of what 
is observed as the coronal source Fig. |2] indicates a precipita- 
tion distance of about 13 arcsec from the site of acceleration to 
the footpoints. This value was used for the computation of the 



References 

Aschwanden, M. J., Brown, J. C, & Kontar, E. P. 2002, Sol. Phys., 210, 383 

Battaglia, M. & Kontar, E. R 2011, ApJ, 735, 42 

Boyer, R., Sotirovsky, P., Machado, M. E., & Rust, D. M. 1985, Sol. Phys., 98, 

255 
Brown, J. C. 1971, Sol. Phys., 18, 489 
Brown, J. C. 1973, Sol. Phys., 31, 143 

Brown, J. C, Aschwanden, M. J., & Kontar, E. R 2002, Sol. Rhys., 210, 373 
Canfield, R. C. 1974, Sol. Rhys., 34, 339 
Chen, Q. R. & Ding, M. D. 2006, ApJ, 641, 1217 
Ding, M. D., Fang, C, & Yun, H. S. 1999, ApJ, 512, 454 
Emslie, A. G. 2003, ApJ, 595, LI 19 
Fang, C. & Ding, M. D. 1995, A&AS, 110, 99 

Fletcher, L., Hannah, I. G., Hudson, H. S., & Metcalf, T. R. 2007, ApJ, 656, 1 187 
Hannah, I. G., Kontar, E. R, & Sirenko, O. K. 2009, ApJ, 707, L45 
Hudson, H. S. 1972, Sol. Rhys., 24, 414 

Hudson, H. S., Wolfson, C. J., & Metcalf, T. R. 2006, Sol. Rhys., 234, 79 
Hurford, G. J., Schmahl, E. J., Schwartz, R. A., et al. 2002, Sol. Rhys., 210, 61 
Kontar, E. R, Brown, J. C, & McArthur, G. K. 2002, Sol. Rhys., 210, 419 
Kontar, E. R, Hannah, I. G., Jeffrey, N. L. S., & Battaglia, M. 2010, ApJ, 717, 

250 
Kontar, E. R, Hannah, I. G., & MacKinnon, A. L. 2008, A&A, 489, L57 
Kontar, E. R & Jeffrey, N. L. S. 2010, A&A, 513, L2+ 
Kretzschmar, M. 2011, A&A, 530, A84+ 

Krucker, S., Hudson, H. S., Jeffrey, N. L. S., et al. 2011, apJ, in press 
Lemen, J. R., Title, A. M., Akin, D. J., et al. 2011, Sol. Rhys., 106 
Lin, R. R, Dennis, B. R., Hurford, G. J., et al. 2002, Sol. Rhys., 210, 3 
Machado, M. E. & Rust, D. M. 1974, Sol. Rhys., 38, 499 
Matthews, S. A., van Driel-Gesztelyi, L., Hudson, H. S., & Nitta, N. V. 2003, 

A&A, 409, 1107 
Meegan, C., Lichti, G., Bhat, R N., et al. 2009, ApJ, 702, 791 
Metcalf, T. R., Alexander, D., Hudson, H. S., & Longcope, D. W. 2003, ApJ, 

595, 483 
Najita, K. & Orrall, R Q. 1970, Sol. Rhys., 15, 176 
Neidig, D. F 1989, Sol. Rhys., 121, 261 

Neidig, D. F, Wiborg, R H., & Gilliam, L. B. 1993, Sol. Rhys., 144, 169 
Potts, H., Hudson, H., Fletcher, L., & Diver, D. 2010, ApJ, 722, 1514 
Saint-Hilaire, R, Krucker, S., & Lin, R. R 2010, ApJ, 721, 1933 
Schmahl, E. J., Pernak, R. L., Hurford, G. J., Lee, J., & Bong, S. 2007, Sol. Rhys., 

240, 241 
Su, Y, Holman, G. D., & Dennis, B. R. 2011, ApJ, 731, 106 
Sylwester, B. & Sylwester, J. 2000, Sol. Rhys., 194, 305 
Svestka, Z. 1970, Sol. Rhys., 13, 471 



M. Battaglia & E. P. Kontar: Multi wavelength source heights 



Vernazza, J. E., Avrett, E. H., & Loeser, R. 1981, ApJS, 45, 635 
Wachter, R., Schou, J., Rabello-Soares, M. C, et al. 2011, Sol. Phys., 19 
Watanabe, K., Krucker, S., Hudson, H., et al. 2010, ApJ, 715, 651 
Xu, Y., Cao, W., Liu, C, et al. 2006, ApJ, 641, 1210 



